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A SUMMARY AND EVALUATION OF SEMI -EMPIRICAL METHODS 
FOR THE PREDICTION OF HELICOPTER ROTOR NOISE 

By 

Robert J. Pegg 
SUMMARY 

A compilation and evaluation of existing semi -empirical methods for 
predicting various helicopter rotor noise components is given in this report. 
These prediction methods should be useful to preliminary design groups, planners 
and user/operators requiring easy-to-use techniques for the evaluation of the 

acoustic characteristics of helicopters. 

The author has selected appropriate available theoretical methods which 
can be solved easily through closed-form solutions or with the help of computer 
analyzed graphs. Examples are included which show the use of the prediction 
techniques and are compared to experimental data. 

INTRODUCTION 

The far-field noise produced by rotors is comprised of periodic compo- 
nents; a series of tones occurring at the blade passage frequency and multiples 
thereof, and broadband components; a random noise arising from inflow from 
inflow turbulence and wake-shedding induced force fluctuations (see figure 1). 
In addition to these primary noise components, rotors can sometimes produce 
an impulsive tvpe noise generally called blade slap which arises from either 


bUde vort«K Interaction or compressibility (thickness) effects or both. The 
complex aerodynamic flow field associated with a helicopter stronoly Influences 
noise components which are Important In different portions of the flight enve- 
lope. During level flight with a low advancing tip Mach number, the rotational 
broadband noise can be considered as the major noise components. At high 
advancing tip Mach numbers, thickness and other compressibility effects can 
be expected to predominate. When the helicopter Is flying In certain flight 
(particularly descent) conditions, blade-vortex Interactions can occur and 
and may result In Intense impulsive noise. 

The purpose of this report is to provide a compilation of existing semi- 
empirical methods for predicting the various rotor noise components for given 
flight operating conditions. These prediction method^ should be useful to 
preliminary design groups for parametric studies and for extrapolating to new 
designs; to planners who want to compare one type of operating procedure with 
another or want to compare the noise signatures of various helicopters for a 
given procedure; and to the user/operator requiring an independent evaluation 
of the acoustic characteristics of several helicopters. 

In the text of the report, the empirical prediction procedures are 
described. The description includes the input and output parameter list, 
steps of computation, required equations and graphs, and the range of validity 
for each part of the prediction procedures. A section is also provided 
showing illustrative examples of the methodology developed in the report. The 
examples include calculations for rotational noise, blade/vortex interaction 
noise, and high-speed thickness noise. 

Appendices, in which considerable attention has been devoted to the 
chronological development of the rotor blade noise technology, have been 
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included l„ effort to provide the user ulth . greeter Insight Into the use 
Of the prediction methods. A diagram of the various retor noise con„onents 
Is Shown in figure 1. The sequence presented In this figure has been followed 
in the general organization of the material In both the text and the Appendices. 
The material In the Appendices has the additional advantage of giving the 
reader some Insight Into the reasoning for reco^mndlng certain prediction 


procedures and the options available to the reader If variations In the 

pr«11ct1on methods must be considered. Two recent publications, reference 1 

and 2, present helicopter noise calculation procedures which are similar to 

the level of prediction sophistication outlined In this report, but which 

are primarily graphically oriented and do not allow the calculation flexibility 
Of this report. 



SYMBOLS 

blade area, 
correlation area, 

blade tip area affected by vortex, 
number of blades 
blade chord, m 

speed of sound at sea level standard conditions, m/sec 
amplitude of first load harmonic 

‘l.T- S.c* S.T loading coefficients In the rotational 

noise equation (A-4) 

incremental blade drag coefficient 


rotor diameter, m 


Da L. rotor disk loading, n/m^ 
mean lift coefficient 
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E 

f 

F 

h 




number of Interactions per revolution 

frequency, Hz 

peak frequency, Hz 

point force 

blade thickness, m 

sunrotion Unit In compressibility noise equetlon 
Bessel function 


thrust constant, N/ro^ 

turbulent boundary layer area correlation constant 
correlation length, m 

blade lift, Newton 

sound pressure harmonic number 

drag divergence Mach number 

effective Mach number = 

1 - Mp cos$ 

forward velocity Mach number 

component of the relative Mach number In observer direction 
blade tip Mach number 

loading fall-off exponent 
rotational speed, rev/sec 
power, watts 


complex sound pressure due to the first A loadlnp harmoni, 
souno pressure due to high freqiency in plane loads. H/n? 

sound pressure due to high frequency drag loads, N/m^ 


1’^ sound pressure due to high frequency thrust loads. H/m^ 

2 

q dynamic head, N/m 

r distance from source to observer, m 

R rotor radius, m 

Rg effective blade radius - xR, m 

RN Reynol ds Number 

^1/3 normal 1 zed broadband noise spectrum 
t time, sec 

T total system thrust, Newton 

u gust magnitude, m/sec 

U frees tream total velocity, m/sec 

V forward velocity, m/sec 

Vj descent velocity, m/sec 

Vj blade tip speed, m/sec 

w rotor induced velocity in forward flight, m/sec 

Wg rotor 'induced velocity in hover, m/sec 

W total acoustic power, watts 

X fraction of radius at which rotor forces are assumed to act 

x^, y^ observer source location coordinates 

a angle between flight path and rotor disk plane, deg. 

$ observer angle measured from the rotor disk plane, deg 

Xg blade loading spectrum function 

it 

6 boundary layer thickness 

AL incremental blade lift due to vortex interactions, Newton 

ASPLbv blade-vortex interaction noise component sound pressure level 


ASPL, 


ASPL. 


ASPL| 

ASPL, 


mr 


ASPL 

Aili 


tr 


Y 


\ 

A 


b> 


'1 


nom» 11 z^ con,p«ss 1 bU 1 ty- 1 nd«ed profll. dPdj nols, c«.^onent 

sound pressure level 

nom.Hzed thickness noise sound pressure level correction factor 

nomellted mein rotor noise conponent sound pressure level 

"or«.lleed thickness noise coniponent sound pressure level 

normollted tell rotor noise component sound pressure level 

•slnkith portion of disk within which b1.de vortex Interaction 
occurs* deg 

enple between flight path and horizontal, deg 

loading harmonic number 

summation limit In equation (A-4) 

sinusoidal gust frequency, rad/sec 

rotor azimuth angle, deg 

complete rotor azimuth of 360*^ 

load solidity 

angle between negative thrust axis and line fro™ hub to observer, deg. 


RECOMMENDED SEMI-EMPIRICAL PREDICTION PROCEDURES 
The background for prediction techniques of helicopter rotor noise 
c^onents Identified In figure 1 has been detailed In the Appendices. Onlv 
e steps necessary for computation are described In this section. These 
prediction techniques are. In some cases. sen,1-emp1r1ca1 which allows their 
ese for the fast calculation of helicopter rotor noise but which may limit 
their universal applicability. So»» degree of engineering insight Is necessor 
to adequately utilize these procedures. Several examples are given In the 
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following section Illustrating the use of the equations and showing the 
correlation with measured noise levels where possible. 

Periodic Noise 

Rotational nolse/steadv and unsteady loading .- The recommended procedure 
for calculating the rotational noise of a helicopter Is the one developed by 
Lowson and Ollerhead (ref. 3). The procedure does not predict the discrete 
noise generated by blade thickness or blade/vortex Interaction. This method 
is based on an arbitrary point loading with random phasing and ..n assumed 
ratio of thrust to drag to radial force of 10:1:1. These force ratios are, 
in reality, varying with forward speed, configuration changes, etc. The 
unsteady load on the rotor Is represented by a spectrum of loading harmonics 
in the computation cf Lowson and Ollerhead. In the high frequency range, 
the loading harmonics are assumed to decay with an exponent of n * -1.85. 

This is a representative value based on the work shown in references 3 and 4. 
The prediction procedure is given as follows. 

Input parameters: 

B number of blades 

m acoustic harmonic number 

Mp forward velocity of helicopter 

blade tip Mach number 

N rotor speed; rev/sec 

r distance from souce to observer 

R rotor radius 

R. effective rotor radius 

e 

T total rotor thrust 
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6 observer angle measured from rotor disk 

Ik rotor azimuth angle 

Output parameters: 

f frequency at a given acoustic harmonic number 

^^^mB sound pressure level at a given harmonic number 

Computational Procedure: 

Step 1. Compute the effective rotor tip Mach number by using: 

+ Mp sinik 

M . -S f (,j 

• 1-MpCos6 

Step 2. Compute the blade passage 'requencles according to: 
f . -JSBN 

1 - Mp cose ' 

step 3. Determine from figure 2 according to the value of 

M^ and 8. The value of the last two parameters must be within the range 
of 

0.5 < Mg < 0.9 ; 0. < B < 80. 

Figure 2 is obtained from reference 3. 

Step 4. Compute the sound pressure level SPLj^jg for each value of 
mB by using: 

SPl^ ■ 20 >09 T ♦ 20 109 -W + m (3) 

pc. 

This procedure can also be used to calculate the rotational noise 
of a tall rotor. Because of the complicated flow environment In which 
the tall rotor normally operates, a different set of Incremental sound 
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pressure levels# ASPL^^* should be determined In Step 3 uslns velues of 
ASPL^^ provided In figure 3. These values are based on an assumed ratio 
of thrust to drag to radial forces of 10:1:0 and a loading exponent 
n ° -2.1. Figure 3 applies for the tall rotor of a hovering helicopter. 
According to reference 5, a high frequency loading exponent n of 
approximately -1.85 can be used for Inflight tall rotor rotational noise. 
Hence# figure 2 may be used for these latter cases. The difference In 
the assumed values of thrust to drag to radial force ratios between 
figures 2 and 3 should be noted In practical applications. 

Rotational nolse/compressi bill tv- Induced profile drag noise .- This noise 
source Is attributed to periodically fluctuating torque loads on a helicopter 
rotor blade generally occurring In high-speed forward flight. The simplified 
prediction procedure presented here Is based on the work of Arndt and Borgman 
(ref. 6). 

Input parameters: 

c rotor blade chord length 

Mg effective blade tip Mach number 

^dd divergence Mach number 

r distance from source to observer 

R rotor radius 

R^ effective rotor radius 

e 

$ observer angle measured from the rotor disk 

Output parameters: 

SPLjj^b component sound pressure level at given blade harmonic 


frequency 


Computational Procedure: 

Step 1. The effective tip Mach number Is obtained according to 

equation (1). If Is less than or equal to "dd- the compressibility- 
induced profile drag noise components vanishes. Proceed with Steps 2 and 

3 If Mg Is greater than Mdd- 

Step 2. Determine ASPL^. from figure 4 by using values of mB. Mg, 
and 8. The range of Mg Is restricted to 

Mdd < Mg < 0.95 . 

Step 3. Compute the compressibility-induced profile drag noise 
component by using the following equation: 

SPl-mB • 20 log -~ + 20 log [(Mg - M^^) c/R] + ASPL^ - 21.6 (4) 

Thickness noise .- The method chosen here for the computation of helicopter 
blade thickness noise was developed by Hawkings and Lowson (ref. 7). 

Input parameters: 

B number of blades 

c blade chord length 

h blade thickness 

Mg effective tip Mach number 

r distance from source to observer 

R rotor radius 

6 observer angle measured from the rotor disk 

Output parameter: 

SPL|nB sound pressure level at a given harmonic frequency. 
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Computational procedure; 

Stepl. [leterml.^e «Sure 5 wcordinj to «1ven values 

of a«, Mj, and d. The range of is confined to 

0.8 < Mg < 0.95 . 

Step 2. Determine ASPL^ from figure 6 according to the hermonic 
number mB and the blade chord to rotor diameter ratio, c/D. 

Step 3. Chfmpute the sound pressure level for the thickness noise 
conponent by using the following equation: 

^09 V 2^" + 20 log h/c + 20 log B + ASPU + ASPL - 

(5) 

In this coinputational procedure, the curves in figure 5 arc coaiputed 
according to the Hawkings and Lowson theory with a value of c/D ■ 0.03. 
The incremental sound pressura level adjustment aSPl^ is a correction 
factor for rotors having different c/0 values. This correction factor 
Is assumed to be independent of observer angle. 

In teraction noise/hlade and vortex interaction. . The computation of nois. 
resulting from the Interaction of a helicopter rotor blade and a trailing tip 
vortex is based on the work of Wright (ref. 8). This noise mechanism Is 
present only for some helicopter flight conditions. Figure 7 illustrates an 
example of helicopter flight conditions where the wake trajectories interact 

with the rotor blades to produce the impulsive noise which is typical of this 
noise component. 

Input parameters: 

B number of blades 

Lq average blade lift 


N rotor speed, rev/sec 

f distance from source to observer 

T total thrust of rotor 

S observer angle measured from rotor disk 

incremental lift on blade owing to Interaction 
total disk azimuth angle, 360“ 

P air density 

Output parameter: 

SPL^B pressure level vt given harmonic frequency 

Computational procedure: 

Step 1. Determine ASPL^^y from figure 8 for g1v<“ v jrfi 

and Arp/ijf^. The range of the aximuth ratio given In f-cure 8 is 

0.02 < < 0.05 . 

Step 2. Compute the sound pressure level for the blade-vortex 
Interaction noise by using the following equation: 

SPL,^ ■ 20 log (r t 20 log -li + ASPU„ . 190.2 (6) 

It has been found from experience that for the noise resulting from a 
blade-vortex Interaction on the advancing side of the rotor disk, a typical 
value of the azimuth ratio Is 0.05. The lift ratio AL/L^ varies 

from about .15 for a conventional blade (constant chord, linear twist, 
square tip) to about .07 for a blade with blade tip air Injection. In 
the present computational procedure, values of the azimuth ratio and the 
Incremental lift ratio are empirically determined. 


Interaction nolse/lnstalldtlon effect .- Increased attention Is being 
directed at subject of noise generated by the Interaction of a turbulent a1r> 
stream with a lifting rotor. There are currently no verified procedures avail- 
able for the calculation of helicopter Installation noise. 

Broadband Noise 

Broadband noise Is a complex noise source comprised of a number of diffe- 
rent components. Current prediction techniques for boundary layer, turbulent 
Inflow, and vortex Interaction effects need additional Investigation to provide 
reasonable correlation with experimental data. The recommended prediction pro- 
cedure Is based on a d1ir«ns1onal analysis approach. Contributions from various 
components, as shown In figure 1, are not explicitly Identified In this 
prediction procedure. 

The recommended rotor-generated broadband noise computation procedure Is 
based on the work of Lowson (ref. 9), Hubbard (ref. 10), Schlegel (ref. 11), 
and Munch (ref. 12). 

Input parameters: 

Ajj total blade area of the rotor 

average blade lift coefficient 

T total thrust of the rotor 

Vj rotor tip speed 

r distance from source to observer 

01 angle between negative thrust axis and line from hub to 

observer, deg. 
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Output parameter: 

SPL ,^3 one-tMrd-oct»ve band sound pressure level for the overall 

contribution from all broadband noise components. 
Computational procedure: 

step 1 . Compute the peak broadband noise frequency by using the 
following equation*: 


fp » -240 log T + 2.448 V-j. + 942 


(7) 


step 2 . Identify the peak frequency with a standard l/3-octave band 

Step 3. compute the peak l/3-octave band sound pressure level by 
using the following equation: 


SPL 


1/3 * 


^0 log (J) ,0 log ^ (cos^ . . 1 ) .s,,, t f(e^) . ,30 


( 8 ) 

where 5,^3 can be determined from figure 9. The lift coefficient 
function. f(C^). for most helicopters with < .48 is 10 log (E^/. 4 ) 
Preliminary experimental data indicate significantly higher values of 

f(CL> for C^ • .48. The normalized spectrum function 5,^3 as shown 
In figure 9 Is developed by Schlegel (ref. 11 ). 

step 4. Compute the sound pressure level for all the other one- 
third-octave bands according to equation ( 8 ) and figure 9 . 


•This equation is dimensionally Incorrect and it is unit dependent, for the 
English system of units this equation should read: 
fp = -240 loq T + .746 Vj + 786 . 
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APPLICATION AND EVALUATION 


This section of the report Illustrates the use of the recommended predic- 
tion methods by numerical examples. In addition, an Indication of parameter 
sensitivity and the effect of variation In any empiricism are Investigated. 

The examples chosen represent a wide range of operating conditions and vari- 
ables. The particular examples show how to apply the prediction techniques 
and the corresponding agreement with existing experimental data. The first 
example specifically Illustrates the calculation of rotational noise from the 
main and tall rotors of a large hovering transport-type helicopter. The 
second example Illustrates the prediction of Impulsive noise resulting froir 
bl ado/vortex Interaction for a wind-tunnel mo(tel In a simulated rate of 
descent. The calculation of high-speed Impulsive noise Is Illustrated In the 
third example and compared with experimental wind-tunnel model data. 

Example 1 

The noise from a large, hovering, single-ma In-rotor helicopter can be 
calculated using equations (2), (3), (7), and (8). Only rotational and broad- 
band noise components are considered for this example because the tip speed 
Is sufficiently low so that thickness and compressibility noise do not affect 
the acoustic spectra. Since the helicopter Is hovering no blade/vortex Inter- 
action Is as'^^d to occur. For the conventional type helicopter, this pro- 
cedure must also be repeated for the tall rotor and the results superimposed 
on a frequency basis. The following Information Is required to predict 
the noise: 
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Main Rotor 

.6 


Tall Rotor 

.6 

21.1 


Tip Mach number 
Rotor speed, rps 
Blade radius, m 
Observer angle, i?, deg 
Number of blades 
Thrust, N 
Blade area, m^ 

Air density, kg/m^ 

Observer distance, m 
^tean lift coefficient 


3.5 

9.45 (31 ft) 
5® 

5 

69420 (15,600 lbs) 
18.58 (200 ft^) 
1 228 

(.00238 slug/ft^) 
61.6 (202 ft) 
.438 


1.52 (5 ft) 
80“ 

5 

5206 (1,170 lbs) 
3.44 (37 ft^) 

1 *228 

(.00238 slug/ft^) 
62.2 (204 ft) 
.182 


The tall rotor rotational noise and thrust Is based on the assumption 
that the ratio of rotor blade thrust to drag Is 10; 1 and that the distance 
from the tall rotor to main rotor Is 11.1 m. The radius and Mach number 



at which the rotor forces are assumed to be concentrated Is .9 of the 
actual radius. 

The rotor noise for the example hovering helicopter is composed of both 
rotational and broadband components for the two rotors. The rotational noise 
Is calculated first from equation (3) and figure 2. Values of the tip Mach 
number are obtained by Interpolation In figure 2. 


Rotational ; » 20 log + 20 log 


69420 


(1.228) (344)^. 45) 


+ ASPL_^ 
\c mr 


= 17.2 - 45.4 + ASPL 

mr 

» -62.6 + ASPL 

mr 

fundamental frequency » s 17.5 hz 

(1 - Mp cose) 


mB 

Frequency 

ASPL„ * 
mr 

Calculated, 

Experimental , 

5 

17.5 

148.4 

85.7 

85.0 

10 

35.0 

143.5 

80.8 

79.5 

15 

52.5 

140.2 

77.5 

76.0 

20 

70.0 

138.4 

75.7 

73.0 

* From figure 3 for Mg 

* .54 and 3 

0 

0 

00 

It 


** Data 

from reference 

13, figure 17. 




17 


Tall rotor; 


20 log ♦ 20 log I2S6 * 4SPI,, 

(1.228)(3M)^(1.52)^ 


• -33.1 - 36.2 + ASPL^^ 

- -69.3 + ASPL^^ 

fundamental frequency * -iiliiLD — s 105.5 Hz; 

l-(0) cos 80 


mB 

frequency 

ASPL^^* 

Calculated, 

Experimental 

5 

:j5 

154.4 

85.1 

82 

10 

210 

149.2 

79.9 

76 

15 

315 

145.9 

76.6 

75 

20 

420 

143.8 

74.5 

73 

*From figure 3 for 

e 

, = .54 and 

It 

00 

o 

o 


**Data 

from reference 

13, figure 17. 




Broadband : 

The rotor broadband noise is calculated based on equation (8). For ' 
the main rotor: 

^^** 1/3 “ 20 log (|JJ) + 10 log (cos^ 85 + ‘^)j + 

+ 10 log + 130. 

= -13.1 - 32.8 + Sy^ + .39 + 130. 

= 84.5 ♦ S^3 
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For the tail rotor: 

SPL^/3 « 20 log (|J|) + 10 log (cos^ 10 ♦ .1)J ♦ $^3 + 

♦ 10 log + 130 

» ”13.9 - 30.2 + S^3 - 3.42 + 130. 

• 82.5 + Sy3 

The pectrum correction value, 5^3* is obtained from figure 9. Peak fre- 
quencies for both main and tail rotors are calculated from equation (7): 
fp « -240 log (69420) + 2.448(208) + 942 

« -1162 + 509 + 942 = 289 Hz 

fp « -240 log (5206.5) + 2.448(202) + 942 
« -892 + 494 + 942 = 544 Hz 

The third octave frequency spectra for the broadband noise uses the above 
information and is calculated below for both the main and tail rotors. The 
main rotor broadband peak frequency of 289 Hz falls within a band where the 
center frequency is 315 Hz, therefore: 


1/3 Octave 
Center Freq. 

^1/3 


40 

-19.5 

64.8 

80 

-15.3 

69.0 

160 

-11.7 

72.6 

315 

- 7.5 

76.8 

630 

-11.5 

72.8 

1,250 

-12.1 

72.2 

2,500 

-16.5 

67.8 

5,000 

-17.0 

67.3 


Tall rotor peak frequency* fp ■ 544 Hz, lies In a third octave band with a 
center frequency of 630 Hz* therefore: 


1/3 Octave 
Center Freq. 

^1/3 

S'’h/3 

20 

-29.0 

56.1 

40 

-24.5 

58.6 

80 

-19.5 

63.6 

160 

-15.3 

67.8 

315 

-11.7 

71.4 

630 

- 7.5 

75.6 

1*250 

-11.5 

71.6 

2*500 

-12.1 

71.0 

5,000 

-16.5 

66.6 

10,000 

-17.0 

66.1 


Discrete noise components at blade passage frequencies previously calculated 
may be combined with the one- third-octave broadband noise levels by adding 
the sound pressure levels in each 1/3-octave band. The calculated values are 
then combined and shown in the following summary table. The table is carried 
out to a frequency where most discrete data for conventional helicopters 
become less than the broadband noise. 
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1 /3-Octave 
Center Freq 

~6 

20 

25 

31.5 

40 

50 

63 

80 

100 

125 

160 

200 

250 

315 

400 

500 

630 

800 

1 K 

1.25 K 

1.6 K 

2 K 
2.5 K 


Tall rotor. SPL 
Rot. B. B. 


76.6 

74.5 


63.6 


Main rotor, spi 
Rot. B. B. 


64.8 


Total, 

SPL 

85.7 

54.1 

80.8 

66.5 

77.5 

75.7 

70.1 

85.1 

73.8 

79.9 


Correlation of the celcelateO soend pressure levels ind those flight test 

»«esur«nents in reference 13 Indicate good agt^e^nt to the fourth hamonic 
for both the main and tail rotors. 

A sensitivity analysis of specific parameters in the rotational and hr 
bend noise components for the previous exan^le that are difficult to measur 
experimentally indicates that a 10 percent variation in: a) thrust will 


result In approximately .9 dB change In noise level and 10 Hz shift In 
broadband frequency; b) air density will change the level approximately .9 
dB; c) speed of sound will change the level approximately 1.8 dB and; d) 
rotor speed will change the broadband level approximately 2.8 dB and shift 
the center frequency 50 Hz. 


Example II 

This example Illustrates the use of the blade-vortex Interaction equation 
(eq. ( 6 )), as applied to a two-bladed rotor wind-tunnel model. The experi- 
mental data used for correlation in this example are reported In reference 14. 
The contribution from the rotational and thickness noise components Is not 
affected by the periodic loading noise Induced by the blade/vortex Interaction 
and Is therefore not considered In this example. The necessary parameters 
required to determine the discrete harmonics resulting from a blade/vortex 
Interaction and the specific values for the example are listed below: 


Total rotor thrust, N 

458.35 

(103 lb) 

Tip Mach number 

.433 


Rotor speed, rps 

22.75 


Observer angle, deg 

490 


Observer distance, m 

1.02 

(3.35 ft) 

Forward speed, m/sec 

21.3 

(70 ft/sec) 

Speed of sound, n^sec 

344. 

(1130 ft/sec) 


It Is assumed that the major contribution of the Interaction occurs over 
that portion of the azimuth which starts when the vortex first Intercepts the 
tip of the blade and the blade sweeps through approximately 18«. The momentary 
blade load disturbance* AL, caused by this Interaction Is assumed to be 


20 percent of the loed of one blade or 10.3 1b. 
hannonlc levels from equation (6) and figure 8 


The rms sound 
a»*e therefore 


pressure 


'®9 (.2 cos 49) * 20 log 1958.351(22.70' 

(1.228)(344)^1.02) * “’’‘‘Ov * 

• -17.6 - 73.8 t dSPt^^ + igo.2 

• 98.8 + 




These calculated harmonics make up an envelope and -re suoerln.. a 
— ntal data taken fr. I - 7 - - 

zt:t " "" -- - — - u 9I . 

The ha JnT V"*'* '-‘Pe - similar. 

Is due to the T passage frequencies 

to the steady .tatlonal and blade thickness sou.es The „ . a. 

““ •’ “ »»-» «. ~ «»., :: ;:i:. 
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Example III 

This example illustrates the prediction and correlation of the high tip- 
speed noise sources. Experimental data for this example are taken from 
reference 15. These data were taken from a hovering 1/7-scale UH-1 rotor 
operating at two thrust conditions. Comparisons of the calculated results 
and the experimental data are given in figure 11 for two tip Mach numbers. 

The following information is required for the calculations: 

B « 2 

r « 3.14 m (10.3 ft) 

6 » 0 * 

T = 485 N. 676 N (109 lb, 152 lb) 

Mg « .9 (N » 47.2 rps, f^ « 94 Hz) 

Mg = .8 (N = 41.95 rps, f^ » 84 Hz) 
c = .0762 m (.25 ft) 

R = 1.05 m (3.43 ft) 

M^jj = .8 (fig. 15) 
h/c = .12 

Rg * .945 m (3.09 ft) 

Numerical calculations will be made to compute compressibility noise, 
rotational noise and finally the contribution of thickness noise. These calcu- 
lations are then compared with the measured data. The compressibility-induced 
profile drag noise is determined from equation (4) and figure 4 for M » .9: 

SPL « 20 log ^ + 20 log [(.9 - .8) .0762/1.05] + ASPU - 21.6 

c 

» -10.4 - 42.8 - 21.6 + ASPL. “ 74.8 + ASPL 
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nd) 

frequency f 
Hz 

ASPlj 

(from eq. (3)) 

2 

94 

173.9 

99.1 

4 

188 

175.3 

100.5 

8 

376 

176 

101.2 

12 

564 

175.7 

101). 9 

16 

752 

174.9 

100.1 

20 

940 

173.9 

99.1 

30 

1410 

170.2 

95.4 


These values are shown in the lower portion of figure 11 as the square data 
points. 

The compressibility-induced drag noise calculated values for M = .8 are 
zero because the tip Mach number is equal to the drag divergence Mach number. 

Thickness noise values are similarly determined from equation (5) and 
figures 5 and 6. 

SPL^ ' 40 log (.9) + 20 log + 20 log (.12) + 20 log (2) + 

+ASPLj + ASPL,^ - .9 
= -24.6 + ASPLy + ASPL,^ (Mg = .9) 

» -26.7 + ASPL^ + ASPL,^ (Mg « .8) 
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mB 

freq. 

Hz 

ASPL|^ 

(for M“.8 & .9) 

4SPLy0Mj«.8 

SPL.dB 

ASPLj0Mj».9 

SPL.dB 

2 

94 

3.5 

117.9 

94.7 

119.6 

98.4 

4 

188 

3.5 

122.1 

99.6 

119.6 

104.1 

8 

386 

3.4 

123.9 

100.6 

129.7 

108.4 

12 

564 

3.2 

123.2 

99.7 

131.6 

110.1 

16 

752 

3.1 

121.7 

98.1 

132.4 

no.8 

20 

940 

3.0 

119.6 

95.9 

132.7 

111.0 

30 

1410 

2.5 

113.4 

89,2 

132.2 

110.0 


The rotor is producing thrust and therefore a rotational noise is also 
calculated. This is done using equation (3) and figure 2 for ».n effective 
tip Mach number of .72 and .81 (Mach number at assumed force location .9R): 



SPhnB- 

8 

.. . (.945) . 

485 


ASPL 


.u 109 „ 228)(344)2(1.05)^ 

-10.4 - 50.4 <■ OSPLj,^ • -60.8 t 

“ nir 

mB 

freq. 

Hz 

ASPL^^@Mg=.81 

(from eg. (2)) 

4SPL„,9H^=.72 

SPL^,dB 
(from eg. (2)) 

2 

94 

155.5 

94.7 

154.3 

93.5 

4 

188 

156.7 

95.9 

154.2 

93.4 

8 

376 

157.3 

96.5 

153.3 

92.5 

12 

564 

156.6 

95.8 

150.6 

89.8 

16 

752 

155.5 

9A.7 

147.7 

86.9 

20 

940 

154.5 

98.7 

145.9 

85.1 

30 

1410 

152.4 

91.6 

142.3 

81.5 
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The sutfinatlon of the three noise component levels for this rotor are 
listed below: 


mB 

frequency, Ht 

SPL 9 Mj » .8 

SPL 0 My » .9 

2 

94 

97.2 

102.6 

4 

188 

100.5 

106.1 

8 

376 

101.2 

109.4 

12 

564 

100.1 

110.7 

16 

752 

98.4 

111.3 

20 

940 

96.3 

111.4 

30 

1410 

89.9 

110.2 


These calculated results are compared with the experimental data of 
reference 15. The calculated acoustic spectrum shapes, shown as envelopes 
in figure 11, are very similar to, the experimental spectrum shapes and the 
individual harmonic levels are generally within 3-8 decibels of the measured 
values. Some of the differences between the calculated and measured values 
in figure 11 can be attributed to the fact that the experimental rotor had 
a NACA 0012 airfoil section (blunt leading edge) whereas the thickness noise 
calculations are made using a circular arc airfoil (sharp leading edge) for 
ease of analysis. 

From figure 11, it can be seen that the primary noise component is that 
from thickness noise. Compressibility and rotational noise fall off more 
rapidly with increasing harmonic number. A ten percent variation in tip Mach 
number will vary the predominant thickness noise harmonic terms approxi- 
mately 1.7 dB. 


CONCLUDING RENARKS 


Semi -empirical methods for the calculation of helicopter main- and tall- 
rotor noise have been presented. These methods represent* In the author's 
opinion, the state of the art for s«n1 -empirical helicopter noise prediction. 
The methods have been presented In such a way that solutions can be obtained 
w1^ ttie help of computer-derived graphs. Examples have been provided 
Illustrating both the analysis procedure and the degree of agreement with 
experimental results. Agreenrent of the analytical and experimental results 
is generally good. The calculation procedures have been organized so that 
as Improvements are made* they can be easily Incorporated In the general 
calculation procedures. 
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Figure 7. Example of flight region where blade/vortex impulsive 
rotor noise may be expected. 
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Figure 9. Normalized rotor broadband noise empirically 
determined spectrum shape. 









theory based on steady loads. 
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Figure 15, Dreg divergence Mech number es e 
function of tip thicknesa for typical 
airfoils used in helicopter rotors. 
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Figure 18. Non-dimensional intense impulsive noise boundary. Level 
flight non-dimensional induced velocity relationship is also 
shov^n. 








APPENDIX A 


PERIODIC NOISE 
Rotational Noise 

The contribution to rotor and propeller discrete noise produced 
by rotating blade forces (thrust, torque, coning) Is called rotational 
noise. This noise Is generated by both the azimuthally uniform steady 
and nonuniform (periodic) forces acting on the blades. Because the 
nonuniform blade forces are harmonic In nature, they are frequently 
referred to as blade-loading harmonics. These forces are the source 
strength of the acoustic dipole portion of the solution to the wave 
equation. Rotation of the blades produces movement successively toward 
and away from the observer, thus Introducing effects due to a varying 
Doppler frequency shift. This modulated Doppler effect causes each 
loading harmonic to produce more than one sound harmonic In the stationary 
system. The level of each sound harmonic Is thereby the result of the 
contribution of many loading harmonics. The periodic repetition of these 
forces at imiltiples of the blade passage frequency results In a discrete 
frequency acoustic spectrum. 

The analytical approach presented in this report Is based on the 
well-known compact source theory. In recent years, more sophisticated 
theories for rotating blade noise from noncompact source distributions 
have been developed. These latter theories require extensive computer 
capability and will not be treated In this report. 
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steady loading theory ,- The earliest attempt to formulate a theory 
of propeller acoustic radiation appears to be that of Lynam and Webb 
(ref. 16)» who recognized that the rotating blades would cause both a 
periodic modulation of the flow through the disk and a corresponding 
acoustic d1stuH>ance at large distances. Lynam and Webb modeled the 
action of the blades by a continuous ring of stationary sources and 
another ring of sinks, with an arbitrary axial displacement; and alter- 
native approach was tried by Bryan (ref. 17) who sought a solution for the 
sound field of a point source In circular motion. Bryan's paper is 
Interesting as an early, though unsuc:essful , example of the retarded-time 
approach. 

Neither of these papers predicts *1 the absolute amplitude for the 
radiated sound, and It was left for Gutin (ref. If) to show that, within 
the frame work of linear acoustics, the steady aerodynamic forces on the 
rotating blades could be represented by acoustic doublets which are a 
function of the torque and thrust of the rotating blades. This develop- 
ment Is limited to a hovering helicopter rotor, and the approximations 
limit the accuracy of the calculations of noise to the fundamental and 
first few harmonics In the far field. Gutin's formulation for the 
rotational noise Is given by equation (A-1) 

K' • S IT - FlJ: I <"* "e <*-’> 

0 6 c 

A comparison between the theoretical results given by this method and 
experimental data for a hovering helicopter rotor (ref. 19) Is shown 
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In figure 12. The comparison In figure 12 Is typical, In that Gutin's 
approach Is usually accurate In predicting the fundamental blade passage 
frequency sound pressure, but underestimates higher hannonic sound 
pressures. The Gutin theory, using steady loads as Inputs, Is useful 
for predicting the low-order pressure harmonics that are most significant 
for structural problems, but Is not adequate for predicting the higher 
order pressure harmonics that are most significant for subjective 
reaction. 

Gutin's work w ^ extended by Hubbard and Regler (ref. 20) to deter- 
mine the amplitude and directivity at points as close as one-chord-width 
from a propeller tip. Several parameters were analytically Investigated, 
and comparisons between theory and experiment were shown to be adequate 
for tip Mach numbers less than one. Garrick and Watkins (ref. 21) 
further extended Gutin's theory to include the effects of forward axial 
speed (vertical climb condition). They synthesized the sound field 
using the pressure field of a concentrated force moving uniformly at 
subsonic speeds. The sound field is given in terms of an integral over 
the rotor disk and, using a simplifying assumption, over an effective 
ring. Each Integral gives expressions for both the near-field and the 
far-fleld sound. The effective ring expression for the far-fleld case 
(eq. (38) in ref. 21) Is: 



mBn 



(A-2) 


Results from equation (A-2) for typical cases show that the direction of 
maximum noise shifts rearward as forward speed Is Increased. 


V 
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Watkins and Durllng (ref. 22) extended the work of Barrlck and 
Watkins to demonstrate the Importance of spanwise and chordwise load 
distribution on the harmonic distribution as well as on the overall 
noise levels. The effect of the chordwise loading distribution Is of 
primary Importance. 

Figure 13 presents a comparison of typical measured and calculated 
rotational noise spectra based on a measured 40 percent chordwlse- 
Integrated, differential blade-loading spectrum and three different 
assumed chordwise-loading distributions (ref. 23). It can be seen that 
all of the distributions provide similar agreement at low frequencies; 
however, as the harmonic number Increases the results begin to diverge. 

In a majority of case studies, the point loading distribution provided 
the 'oest agreement between calculation and measurement. The rectangular 
distribution always provided poor agreement while the sawtooth distribution 
(which should be a more realistic distribution) provided better agreement. 
The exact shape of the chordwise-loading distribution Is particularly 
important at frequencies where compact source assumptions begin to 
break down. However, In most of the cases studied, this assumption was 
not violated and the point loading distribution was adequate. The gener- 
ally good agreement between the measured and calculated rotational noise, 
based on loading measured at eighty percent span, suggests that the 
aerodynamic loading Is representative of the entire rotor. 

Unsteady loading theory .- A technique for obtaining solutions to the 
rotational noise problem when the blade loads are unsteady was advanced 
by Schlegel, King, and Mull In reference 11. This technique resulted in a 
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rr — --- 

‘ » «.. . recun,.,.. puU. ^ 

2'- «• wm .p.H„„u, ^U . 1 

«. * -..lysis «,„ wm «,.t f Gutm. Th. basic M.,t.t1«, to this 
.. -t -. IS tbst biGbec b,™„„tc 1O.01O, Gst. („ > 3,, „.t .v.lUb,. 

11.r solution «, pr.sent«l by Sadler and Lo«v (ref. 24,. 

Tbe rotational noise analysis of Hrigbt (ref. 2S) yields a far-fleld 

7V7” -- - -d „ botb steady al flulT; 

t^^t ^toroue values In a manner sleliar to that of reference 11. 

.<!. oadlng ,s ass»ad concentrated at one redial position. Tbis 

..^ysls provides excellent Insight Into tbe eodal developr^nt of rota- 

noise and a good basic understanding of tbe secbanis. of noise 
generation due to a rotational pressure field. 

Tbe nore widely used prediction tecbn1,ue for rotational noise Is 

0? r — -- 

point force In arbitrary notion was given by: 


«1 - h 

^Cor(l-M^)2 


•• ” 1-M^ « J 


(b-3) 


"h.™ F, represents tbe conponents of tbe point fore., and H Is tbe 

compon6nt of the relatii/A ^ i. ^ ^ 

Tbe relative Hacb number In tbe direction of tbe observer. 

Observer and source locations are *•*#< 

lU. ■ s , .Tionsare .. and y, . respectively, and terns 

was d • ‘•‘^'<1-. 

developed for calculating tbe roUtlonal noise, using tbe rotor 


thrust, the drag and radial components of blade loading, and a compact 
source model. The effect of unsteady loads and the hub convection Hach 
number transformation was Included also. In this method, rotor loading 
harmonics are assumed to decay as Inverse functions (X"^) of the 
loading harmonic number X. The exponent n is determined empirically. 
In addition, the phase between the loading harmonics has been randomized 
to simplify the calculation procedures, that Is, loading harmonic 
pressures are added In the mean-squared sense. 

Equation (A-3) Is the basis for the rotational noise equations which 
are used In the Recoiwnended Prediction Procedures section of this report. 
The equations used here have been derived In reference 1 as described 
below. 

The complex sound pressure due to the first A loading harmonics: 
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The mean-squared pressures due to the higher loading harmonics are: 
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Equations (A-4a) through )A-4d) are then added in a mean-squared sense 


"* X<A "*r X>A '"d X'^A "t X>A 


(A-5) 


The unsteady forces acting on a rotor blade come from many mechanisms; 
rotor/rotor wake interaction, atmospheric turbulence, blade motion, and 
airfoil characteristics, among others. An analytical method that can 
accurately predict the combined effect of all the sources mentioned above 
is still beyond the present state-of-the-art. 


In their studies. Lowson and Ollerhead (refs. 3. 27) noted that 
measured blade airloading harmonic amplitudes tended to fit an exponential 
decay relation with harmonic number. This takes the form: 

“ ^0 ^ ” (A-6) 

where is the amplitude of the first harmonic, X is the loading 
harmonic number, and n is the decay exponent. The loading decay 
rate for a given configuration will vary with the type of flight and 
atmospheric conditions. Increased high frequency airloads can be expected 
to have more effect on the higher order sound harmonics than the lower 
order ones. Unfortunately, very little experimental data are available 
beyond the 10th airloading harmonic so that the airloading harmonic 
spectrum is far from being completely defined. Reference 2 gives a review 
of existing data on loading decay rates for a variety of operating 
conditions and rotors. 

Lowson and Ollerhead (ref. 3) have shown that loading harmonics of 
order X are significant radiators into the m^^ sound harmonic in the 
range: 


mB (1-M) < X < mB (1+M) (A-7) 

It has been determined by Stemfeld et al (ref. 28) that phase 
differences between the lower loading harmonics were of importance. Since 
precise harmonic phasing of the airloads cannot be determined (being to 
a large extent a function of random factors), the logical approximation 
is to assume that each loading harmonic is randomly phased with respect 
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to the others, since this assiaptlon Is not necessarily true for e 
specific situation, relative phasing error betNeen loading hanaonics 
could account for the randon phasing, Lowson (ref. 3> adjusted the 
loading exponent by adding .5. Hosier's data («f. 23) substantiates 
randou phasing In a band t 40 degrees ab«,t 0, with the exception of the 
very lowest hanaonics where more precise phase Infonaatlon seeais Important. 

Co apresslbllltv-induced profile Hr., p,rtodlc blade loading 

can have a gross effect on the radiated noise of a rotating blade system 
>d»n the blade speeds are transonic. The sources of these unsteady blade 
loadings for subsonic blade speeds were discussed In the previous section. 
This section deals with the periodic blade torque-force fluctuation due to 
the compressibility associated with high tip speeds. Noise generated by 
these periodic blade loadings have a tendency to dmslnate the noise signa- 
ture of a helicopter when they occur. This Is cosmonly called Impulsive 
noise. 

The rotor Impulsive noise signature can be characterized by at least 

six vaHables which contribute to the overall sharp-slapping or cracking 

effect. These variables are shown In figure 14. Of the variables shown. 

crest factor and rise time are considered to be of first order Importance. 

Psychoacoustic Investigations currently In progress estimate that the 

crest factor (ratio of peak pressure to rms pressure) must be greater than 

15 dB for the signature to be considered Impulsive. Additional Information 

relating to the subjective aspects of helicopter Impulsive noise can be 

fcund In references 9, 30, and 31 respectively. References 32, 33, and 

34 address the problem of providing a rating for helicopter Impulsive 
noise. 
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Ttar. U , ,„cr„,. 1 „ tta „.u. .f the rotor 

s'Llu 'T'"* co„,res- 

bn ty effects such es the .ppesrsoce of oscm.t.ng shod, wsves oo 

the bUdes e,s, become 1.„ort.nt In this tip speed rsnge, this Is csUed 

co^ressIblHty noise, from the scoustic point of view, the thickness 

»d loading noise „so show a similar behavior In both the Incroase of 

"Oise level and the directivity pattern. However, recent studies by 

Box^n. Schmitz add fu and also by Hanson and fink have shown that the 

non inear effects, resulting from the quadrapole source tem m the 

governing acoustic equation of ffowcs Williams and Hawkings. m,y be 

•ISO important in the calculation of cos^rosslblllty noise 

The highly directional characteristic blade lamulslve 

This Th“ 

and has 1 T' "™ber. 

en entifled experimentally by Cox (ref. 31 ) in full-scale 
L"lc'aT"t*h ’"«*‘'B*t1on (ref. 35) 

Sid T ! “ ^ 

olde Of the rotor disk extending from 75 percent radius to the blade tip 
wos experiencing some degree of shock phenonena. Thus, while there Is 
an aconstlc source effect at high tip speeds which Is partially responsible 

TO h n Of impulsive noise, the advancing blade tip Hach nmrber. 

Which is a function of aircraft fn»M..a«r 4 ... ^ 

^“Tward speed, accounts for the highly 

direotlonal character of this type of noise. Arndt and Borgman (ref. 5 ) 

•*e presented an expression showing the effect of c»press1b111ty drag 
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on the Impulsive noise as: 


a 

mB 


mB ° 
4ir^^ 


M R c ^ 2 /, 1 » 

ig r'>'o Z- 6j sine) 


(A-8) 


j»-«0 


Besides aircraft configuration parameters, this expression Is also 
a function of Incremental drag coefficient. . and the azimuth angle. 
A«. over which the Mach number is greater than the drag divergence Mach 
number. M^^* Figure 15 shows the variation of for typical helicopter 
airfoil sections. Equation (A-8) was derived from the basic rotational 
noise equations with only the drag term being considered. 


Thickness Noise 

Blade thickness noise Is due to the volume displacement of the air 
surrounding the blade, as It moves through the air. This noise Is 
strongly directional In the plane of rotation. Observed sound radiation 
from helicopters In high speed forward flight suggests that as the blade 
tip advances a strong pressure pulse Is generated. The noise from this 
source Increases rapidly with advancing tip Mach number and becomes 
considerably stronger than the rotational noise components. 

Emsthausen (refs. 37. 38) using primarily an experimental analysis, 
studied the characteristics of thickness noise. He emphasized the pre- 
dominance of higher harmonics at blade tip Mach number near 1.0. Deming 
(ref. 38) performed experimental and theoretical work with symmetrical 
blade sections operating at zero angle of attack and developed the first 
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theory for thickness noise. The development of this theory follows much 
the same procedures as Gutin's rotational noise work, however, his 
analysis Is essential for noncompact sources. The airfoil Is replaced 
by a distribution of acoustic monopole sources operating next to the 
Infinite wall. The strength of these sources Is the normal velocity 
distribution on the airfoil. Due to the lack of high-speed computers 
Deming Introduced approximations so that the calculations of the fundamental 
and the first few harmonics are realistic. His calculations are limited 
to static propellers with the observer in the far field. The Diprose 
model (ref. 39) for thickness noise considered distributed radial and 
chordwise blade sources. This modal Is a slight Improvement to the com- 
pact source theory. Based on the original work of Billing (ref. 40), 

Amoldl (ref. 41) formulated a procedure which includes the near- and 
fdr-fl6ld thicknsss noise for compect sources » 

Hore recently, Hawkings and Lowson (ref. 7) presented a theoretical 
analysis of supersonic rotor noise. This analysis results In a noncompact 
solution to the Llghthlll theory of aerodynamic noise. The solution Is 
relatively tractable once the blade configuration Is known. The following 
equation Is valid for hovering conditions: 
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The spectrum shape of the thickness noise component Is directly 
affected by the Mach number and the chord- to-dlameter (slenderness) 
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ratio. The spectnn shape Is detererfhed by the Inteprel of equetton A-9. 
F 19 U« 16 shoes how the verletlon of these two perameters affects the 
shape for Inplane «asurenents. Increasing the slenderness ratio (c/0) 
of the rotor blade does not change the spectnie shape - only the level. 
Changing the Mach number, however, changes both the shape and sound 
pressure level. The results of this equation provide reasonable agree- 
ment with the more sophisticated theory of Farassat (ref. 42) out to 
approximately 40 haraonics. Equation (A-9) is used In the Recoemmnded 
Prediction Procedures section of this report. 

Interaction Noise 

In this section, the generation of discrete noise due to the 
impingement or Interaction of a trailed line vortex on a following lifting 
surface will be discussed. Two cases are of considerable Importance to 
aircraft using rotating blade systems: first, the vortex Interaction 
with either a following blade of a helicopter main or tall rotor which 
generally results In an Impulsive noise condition (commonly referred to 
as blade slap) occurring over selected areas of the flight envelope 

and second, the Installation effects wherein tall rotor Inflow or wake 
are distorted by aircraft geometry. 

B _ 1ade/vortex Interactlorr ,- The generation of lift by a helicopter 
rotor blade results In a complex aerodynamic flow-field. This flow-field 
Is composed of turbulent air and shed tip vortices. Rapid blade lift 
fluctuations generated by the Interact: ,n of a vortex from another blade 
will result In Impulsive noise. Noises of this type, when they occur, 
can be predominant. An accurate analysis of this noise Is dependent on a 
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knowledge of the geometiy and dynamics of the concentrated shed rotor 
wake and on the ability to predict the time-variant blade lift distribution. 

Two principal types of helicopter Interaction noise occur; these 
are due to vortex-impingement with the following blade and vortex-lmplnge- 
ment with the tall rotor. Figure 17 shows, schematically, two orientations 
of a vortex Interacting with a lifting blade. Figure 17 (a) shows a 
vortex with axis parallel to the blade, a situation which represents the 
most severe case. In this case, the pressure distribution along a large 
portion of the blade span will undergo a rapid angle-of-attack change 
which results In an Impulsive noise. In figure 17(b). the axis of the 
vortex core Is perpendicular to the blade-span and only a small spanwise 
section of blade Is affected but over a relatively long period of time. 

The Incremental loading experienced by the blade Is caused by the Induced 
velocity field of the vortex which forces a change In the blade section 
angle-of-attack. The magnitude of this angle-of-attack change is directly 
proportional to the strength of the vortex and Inversely proportional to 
the spacing between blade and vortex. 

The Influence of the shed main-rotor wake on the noise generated by 
the tall rotor may be significant In some cases. The primary operating 
parameters associated with this problem are: 

(a) the advance ratio, y. (the main rotor wake Intersects 
the tall rotor at different azimuth locations because of wake 
sweep angle). 

(b) the direction of tall rotor rotation (under some 
operating conditions the wake vortex can either Increase or 
decrease the relative tall rotor tip speed). 
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(c) the tell rotor tip speed, end 

(d) the Rialfi rotor thrust end essocleted tip vortex strength 
(the rotor loeding coefficient provides an Indication of how fast 
the wake moves vertically away from the main rotor and consequently 
where It Interacts with the tall rotor). 

Accurate calculation of a helicopter vortex wake trajectory Is a very 
complex coiq)utat1ona1 process which can be done only with the aid of a 
computer. The primary reason for this is the interaction of the vortex 
filaments with each other. Two such wake programs are currently available 
and are outlined In references 43 and 44. A less accurate, but still 
useful, representation of the early stage of a helicopter vortex wake 
trajectory evolution can be found In reference 45. 

In an effort to avoid the previously outlined complications, a 
technique has been derived to locate the flight region of maximum impulsive 
noise based on currently available experimental noise data (ref. 46). 
Figure 18 presents the maximum Impulsive noise boundary as a function of 
two non-dimensional velocities. 

The nondimensional velocity variables employed in this presentation 
are 

V^ina + w X (abscissa) 

0 "o 

To convert this nondimensional region to the more familiar rate of 
descent against airspeed presentations, the flight conditions and rotor 
loadings must be known. The following equations can be used to perforii 
the conversion: 
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good eorrolattoi. betMon oxper(»«nt.l and oalculatod data was shown In 
this wort. Fllotas (raf. 51) darivas a sound prassura axprasslon by 
modallng a fintta aspact ratio wing flying at a unifonn spaad o«r a carpat 
of aqually spaced. Infinitely long lines of vortices. Using the method 
developed by Wright (raf. 8). an aguatlon for the har«n1c sound prassura 
level has bean derived. This aquation Is dependent on a Fourier repre- 
sentation of a sinusoidal blade loading occurring once per revolution. 

The prediction method used In this report 1$ based on the following 
formulation by Wright: 

•>mB ■ (c; X, 

"here E t"* "umber of Interactions occurring every revolution. 

Is the load solidity or the fraction of the effective disk annulus 
that the unsteady loading region occupies, is a thrust constant and 
Xj is the blade loading spectrum function. 

A research program using small rotors In a static operating situa- 
tion (ref. 52) revealed the effect of configuration on the shed vortex 
characteristics. It was found during this Investigation that an Increase 
In the number of blades reduced the axial spacing between vortices and 
Increased the angular rate at which they revolve about their centroids 

of vorticity, which caused them to become more unstible and to diffuse 
sooner. 

The effect of collective pitch on vortex strength Is that axial 
spacing between vortices decreases along with vortex strength. The tip 
vortices become Increasiggly unstable and diffuse more rapidly, it was 


observed (ref. 52) that the effect of tip speed on vortex trajectory was 
generally not noticed but that unstable sinusoidal fluctuations occur 
close to the rotor once a strong shock wave forms on the surface of the 
airfoil. This Instability was noticed at Mach numbers above .75 for a 
NACA 0012 airfoil. 

Installation effects .- Rotor noise due to Installation effects 
generally results from the Interaction of the rotor wake Interaction 
with fuselage* wing* nacelle, tall fin, or another rotor. The basic 
mechanisms that alter the rotational noise due to distorted Inflow have 
been previously discussed In the Unsteady Loading Theory section. In 
the case of pusher tall rotor configuration (rotor thrust vector toward 
fin* Inflow over fin) the turbulent wakes from the tall fin can give 
rise to a severely distorted Inflow which significantly Increases the 
noise generated. This occurs 1n a manner similar to the rotor-stator 
Interaction phenomenon (ref. 53) which exists in turbofan engines. This 
unsteady fin force has been measured on a full-scale helicopter In 
reference 54, The effect of interaction of the main rotor wake on the 
blade loading of a tall rotor was Investigated 1n model studies In 
reference 55, It was noted In this experiment that the main rotor wake 
Impingement location on the tall rotor produced Increased tone noise at 
tall rotor blade passage frequency and produced the occurrence of combi- 
nation tones, I.e., tones which are a function of both the main and tall 
rotor rotational frequency. Figure T9 presents acoustic spectra to 
5000 Hz and shows the effect of wake location on the tail rotor noise. 
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BROADBAND NOISE 

IiHtlil 1nvtst1«at1m Into the basic mechanlsns of Isolated airfoil 
noise nas gelded by Sharland (ref. 56) and Curie (ref. 57). He suggested 
that broadband noise fro« a surface In a noving flow can be described by 
at least three mechanlsns which cause random fluctuating forces. One 
fom of random fluctuating force Is the surface pressure field arising 
from a turbulent boundary layer. Another form Is the force fluctuations 
resulting from shed vorticity. Finally, If the surface Is moving In a 
flow which Is Initially turbulent, a third component of fluctuating f-rce 
can arise. Once the mechanisms of noise generation am postulated In a 
physical sense, analytical models can be derived to Indicate how much 
acoustic power Is generated from each source mechanism and which parameters 
affect the source strength. A detailed discussion of these mechanisms 
Is contained In reference 58 and Is suinarlaed next. The relative Is^r- 
tance of the various broadband noise sources are then presented In 
sunmary figure for a typical helicopter rotor blade. The sound pressures 
were calculated along the rotating blade at 10 equal blade areas. 

Inflow Turbulence Noise 

This section will outline the radiation from a plate In a turbulent 
airflow. If the speed Is subsonic, the acoustic wavelength will be 
larger than the turbulence scale so that sound from plates of dimensions 
up to several times the turbulent scale should be calculated simply 
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from the fluctuating pressures. The expression for the fluctuating 
density from a distribution of fluctuating sources can be written as 
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The total acoustic power, W. Is the integral of 
which Is 
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The r.tanM source time correlation of the rate-of-chanje of pressure 
on the plate must nou be evaluated. To obtain equation (B-7) It has 
been ass«« that the plate dimension uas of the orter of the acoustic 
wavelength. If it is assumed that this dimension Is small compai«d to 
the wavelength then the "compacf source assumption can be utlllted and 
the reUrded time assumption dropped. One of the Integrations over 
the surface Is simply a correlation Integral which can be written; 
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»diere A^ Is the correlation area of the pres. ires. Substituting 
equation (B-3) Into equation (B-2) gives: 


w * 


12npc, 


at 


A^dA^ 


(B-4) 


If the blade Is In a homogeneous path of turbulence* then the Integral 
simply becomes a multiple of the blade area* A|^. To evaluate 
equation (B>4) for the case of Impinging turbuience we may take advan* 
tage of the results of Kemp and Sears (ref. 59) as simplified by 
Lowson (ref. 58). This shows that the mean lift per unit span* L* on 
an airfoil In a sinusoidal gust of frequency, u* and magnitude* u* Is 
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where c. Is the airfoil chord. To extend this to three dimensional 
conditions we note that the mean lift In terms of pressure fluctuations 
Is 
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where t Is the correlation length In the chordwise direction, 
c 

Let 
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where a typical fre<;uency, then equation (B*4) becomes 
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The 2A arises because both sides of the plate radiate. Put 

tt> ■ - (B«*9) 

t 

where e Is a typical length and let A^ • » 8 and A » Aj^* 
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the flfiil result for the sound power radiated from an airfoil due to 
Incident turbulence Is 




(B-10) 
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Equation (B-11) has been plotted In figure 19 as a function of Mach 
number for a typical helicopter blade. Brown and Ollerhead (ref. 60) 
provide an approximation to the maximum wake turbulence In terms of the 
lifting surface parameters: 
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This result Is similar to that of Sharland (ref. 56) for the same case. 
The typical frequency for turbulence- induced noise will vary as the 
frequency content of the turbulence that is interacting with the surfaces. 

Boundary Layer Noise 

The same general formula» equation (B-5), can be applied to radiation 
from the boundary layer pressure fluctuations acting on the plate. 

Although theory Indicates that the direct quadrupole radiation from 
the turbulence Is stronger for this case. It is of Interest to attempt 
a boundary layer pressure fluctuation calculations. Equation (B-4) can 
be written as 


B-4 




EuptHnnts indicate that <s «f the order of 36 x 10'* q*. hhert q 

it the dyiUBlc head (J pu*). Hoxewr. there is leas aqreeaant over the 

■agnltudea of the correlation date. Valves of correlation anas assvnad 
by various authors are given In Table I. 


MWUC 1 


Source 


.V 


6 * 


Doak 

144,006*^ 

0.63 U/b) 

58.000 \? 

Lighthlll 

4.406*^ 

0.05 U/(d 

0.011 

Bull and Willis 

1.786*^ 



Shari and 



0.500 

Bull 

1.506*^ 

0.30 U/oi 

0.135 


Most recent studies of boundary layer pressure fluctuations have broadly 
aqreed .ith Bull's data. (ref. 61). Therefore, the correlation area 
can be given by the expression 


where the constant K, is in the range 0.1 < K, < l.o. Substitu ting 
equation (B*14) into (B-13) yields the approximate result. 
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or the sound pressure; 
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Equetlon (B-16) Is shown In'.flqure 20 In comparison with the other 
broadband noise source components. The effect of boundary layer noise Is 
of second order Importance when compared to the other sources. Equation 
(B-15) Is In agreement with Sharland (ref. 56). 


Vortex Noise 

The vortex noise source cen he cetejorlzed into two coovonents. The 

first component results from lift fluctuations produced at the trailing 

edge of the blade by shed vortices, and the second component occurs when 

a tip vortex interacts with a blade tip. An alternate appraach to the 

determination of vortex noise has been presented by Lowson in reference 9. 

This approach is based on a dimensional dependence of sound from fluctuating 
forces . 

Presented below are the relationships for vortex shedding noise as 

hypothesized by Sharland (ref. 56) and tip vortex self-interaction noise 
as stated by Lowson (ref. 58). 

lortgx shedding - Beginning with equation (B-4), an expression can 
be derived for the total radiated power due to vortex shedding. If we 
think of the fluctuating pressure difference as a local lift fluctuation 



pw unit aru. then p can be expressed In terns of the local lift 
coefficient 
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W»re u is the local eean .elocity parallel to the surface. Therefore. 
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in order to describe the Integral In equation (B-ig> „ tenns of the 
flou parameters the mechanisms by uhlch lift fluctuations are produced 
In the absence of external excitation must be knoun. There Is still 
some uncertainty about the magnitude of lift fluctuations Induced by the 
bound circulation changes which accompany trailing edge vortex shedding, 
but It has been suggested that the lift fluctuations can be related to 
time fluctuations In the boundary layer thickness at the trailing edge. 
On this basis, for conditions typical of propeller blades, the root mean 
square of the fluctuating lift coefficient should be of the numerical 
order of the -1/5 power of the Reynolds number. The frequency of the 
lift fluctuations should be about the same as the characteristic fre- 
quency of vortex shedding at the trailing edge: 
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where c Is the chord. For Instance, a rotating blade with a 10-percent 
thickness ratio, the correlation area should be governed by the size of 
the larger eddies at the trailing edge, which right be on the order of 

half the section thickness, that Is .05c^. Equation (B-19) can be 
rewritten as 


48ir c. 


U^o5C^^ AcdA 


or Integrating. 


. . ^4^ . ^ pc V, 


120 c. 


(B-21) 


(B-22) 


A typical sound pressure level assuming a dipole distribution Is 
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J j P . radiatio n.- Utilizing the previously derived equation for 
turbulence Induced noise (eq. B-10)) and assuming that correlation lengths 
are on the order of 5 percent of the chord, a typical frequency Is 


then substituting Into equation (B-17) 
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where A^. Is the blade tip area over which the vortex Interaction is 
to occur. This area is assumed to be approximately 5 percent of the 
blade area. If a dipole distribution is assuned, then the sound pres- 
sure is 
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Recent work presented by Clark (ref, 62) indicates that when blade tip 
twist is used on a helicopter rotor a pronounced reduction in noise level 
is evident. One effect of tip twist is to modify the blade tip vortex. 

The four broadband noise sources; turbulent (eq. B-11)), a vortex 
(eq. (B-23)), tip radiation (eq. (B-26)), and boundary layer (eq. (B-16)), 
are compared in figure 20. This figure shows the calculated variation 
of these noise sources with rotor tip Mach number and rotor radius. From 
figure 20(a), the predominant noise source is that resulting from 
turbulent inflow. The shed vortex noise source is the second source of 
importance. The broadband noise from the boundary layer and vortex 
impingement on the tip are insignificant. Figure 20(b) shows the varia- 
tion of the various broadband noise sources along the rotor radius. The 
primary broadband noise radiation portion of the rotor blade seems to 
be the outboard 50 percent. 

Dimensional dependence of sound from fluctuating forces .- Early 
work on broadband noise performed by Yudin (ref. 63) and later modified 
by Hubbard (ref. 10) and Schlegel (ref. 11) assumed explicitly that the 
sound generation was due to trailing edge vortex effects. This work was 
based primarily on dimensional arguments. The method developed in this 
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section has been correlated with measured broadband noise from helicopter 
rotor blades and Is In reasonable agreement. For a point fluctuating 

force the basic law governing the sound radiation Is well known (Curie, 
ref. 57 ) as 

M (B-27) 

To calcul.ate the sound radiated by a random pressure field acting over 
an area. It Is necessary to split up the field Into correlation areas 
Ag over which the pressure may be regarded as coherent. Thus, the 
equivalent point force Is given by F » A^.P. Because the force Is 
random the overall sound pressure Is given by the sum of the square of 
the effects of the random forces. Thus, If there are E Individual 
events over the blade the total mean square sound pressure Is given by 

(SP)' ' E a / (||)2 


It Is usual to assume that the magnitude of the fluctuating pressures 

acting on the blade Is simply proportional to the mean pressure, that 

Is to the mean thrust divided by the blade area A^. 
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where f Is a typical frequency for the fluctuating forces, in turn. 

f may be put proportional to a typical velocity divided by a typical 
length, t, thus giving 



SI 


The d1«ns1on.l dependence «, hes been used In se*erel dete 

conepses. end W the hesls for the recan«nded p^dtctlon p^cedu^ f.r 

helicopter rotor broadband In this report. The spectral features 

Of the broadband noise are distinctly similar In all cases. It would appear 

to define a “universal" spectral distribution for the broadband noise and to 

relate the center frequency to the aircraft's operating parameters. This 
has been done oy Munch (ref. 12). 

The sound pressure level of a third octave band of rotor broadband 
noise Is 
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et»re 10 logCcos^e, t . 1 ) ,s the directionality effect as derived by lowso 
(ref. 9). The generalized empirical spectrum shape number, s,,,. is shown 
In figure 20. Data collapses were performed using measured rotor data to 
provide this spectrum shape. The peak frequency was determined to be 


fp » -240 log T + 2.448 + 942 

■ -240 log T + .746 + 786 
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The effect of mean rotor lift coefficient has been based on the Initial 
derive ion of equation ,0-31, and empirical flight infp^eion. Broadband 

T h" Hn T" "" 

•«. The lift coefficient function. f(C ). u defined below as- 
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for > .48 


Sound pressure levels for rotors with mean lift coefficients much greater 
than .6 have not been rigorously ve»*1fled. Practical helicopter flight 
conditions, however, generally have lift coefficients less than .6. 
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